The invariants of the local Clifford group 
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We study the algebra of complex polynomials which remain invariant under the action of the local 
Clifford group under conjugation. Within this algebra, we consider the linear spaces of homogeneous 
polynomials degree by degree and construct bases for these vector spaces for each degree, thereby 
obtaining a generating set of polynomial invariants. Our approach is based on the description of 
Clifford operators in terms of linear operations over GF(2). Such a study of polynomial invariants 
of the local Clifford group is mainly of importance in quantum coding theory, in particular in the 
classification of binary quantum codes. Some applications in entanglement theory and quantum 
computing are briefly discussed as well. 



I. INTRODUCTION 

The (local) Clifford group plays an important role in 
numerous theoretical investigations, as well as applica- 
tions, in quantum information theory, quantum comput- 
ing and quantum error correction 2] 3] 4] 5] |fj [7J . The 
Clifford group C\ on one qubit consists of all 2 x 2 unitary 
operators which map the Pauli group Q\ =< a±, ct 2 , c 3 > 
to itself under conjugation, where &\, 02, 03 are the Pauli 
matrices. In other words, C± is the normalizer of Q\ in 
the unitary group U{2). The local Clifford group C l n on 
71 qubits, which is our topic of interest in the following, 
is the n-fold tensor product of C\ with itself. 

In this paper we study the invariant algebra of the local 
Clifford group, defined as follows: let {pij} be a set of 2 2n 
variables, which are assembled in a 2" x 2™ matrix p — 
(pij ) . The invariant algebra of C l n then consists of all com- 
plex polynomials F(p) = F(pu, p\2, . ■ . , P2™2™) which re- 
main invariant under the substitutions p — > U pU^ , for 
every U € C l n It is our goal to construct a generat- 
ing set of this algebra. 

This research started out as the ground work for the 
study of equivalence classes of binary quantum stabilizer 
codes, the latter being a large and extensively studied 
class of quantum codes @ . A stabilizer code is a joint 
eigenspace of a set of commuting observables in the Pauli 
group on n qubits and is described by the projector ps 
on this subspace. Two stabilizer codes ps and ps* on n 
qubits are called equivalent if there exists a local unitary 
operator U € U(2)® n such that UpsU^ is equal to p s > 
modulo a permutation of the n qubits. A natural ques- 
tion to ask is how the equivalence class of a code can be 
characterized by a minimal set of invariants, i.e., (poly- 
nomial) functions F(p$) in the entries of the matrix ps 
which take on equal values for equivalent codes. This 
is, however, a difficult and unsolved problem. Therefore, 
given the explicit connections between stabilizer codes, 
the Pauli group and the Clifford group, it seems natural 
to consider a restricted version of this equivalence rela- 
tion, where only operators U £ C l n are considered, and 
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this is where the invariant algebra of C l n comes into play. 
What is more, it is to date unclear whether this restric- 
tion is in fact a restriction at all: indeed, the question 
exists whether every two equivalent stabilizer codes are 
also equivalent in this second, restricted sense. A possible 
way towards solving this problem is through a study of 
invariants (cfr. also Q)- Moreover, the problem of recog- 
nizing local unitary and/or local Clifford equivalence of 
certain classes of multipartite pure quantum states (sta- 
bilizer states, graph states) has recently gained attention 
both in entanglement theory |3fl [5j :6J and in the one-way 
quantum computing model |lCj . These examples make 
for a number of application domains of the present work. 

From a somewhat different perspective, the invariant 
theory of the Clifford group is also of interest from a 
purely mathemati cal p oint of view. Runge IlJandNebe, 
Rains and Sloane [12| |I3( published a series of papers in 
which they investigate the connection between the invari- 
ants of the (entire) Clifford group (and generalizations 
thereof) and the so-called generalized weight polynomi- 
als of a class of self-dual classical binary codes. Their 
work is a considerable generalization of a central result in 
classical coding theory, known as Gleason's theorem , 
which states that the invariant algebra of C\ is generated 
by the weight enumerators of the class of doubly-even 
self-dual classical codes (the definition of the invariant 
algebra of C\ is here somewhat different than ours, cfr. 
footnote 1). It is interesting that the Clifford group - a 
group which appears naturally in a quantum theoretical 
setting, has such a connection, through invariant the- 
ory, with the theory of classical codes. It is not known 
whether this link is a mere coincidence or a manifesta- 
tion of some deeper result |15| . This remark may serve 
as another justification of the present research. 

In our study of the invariant algebra of C„, we will 
make extensive use of the description of this group in 
terms of binary linear algebra, i.e., algebra over the field 
GF(2) — ¥2- It is indeed well known that n-qubit (lo- 
cal) Clifford operations can be represented elegantly by 
a certain class of 2n x 2n linear operators over F2 [jj |3 
and this binary picture makes the (local) Clifford group 
particularly manageable in the following. In order to ob- 
tain a generating set of the invariant algebra, we will 
adopt the following basic strategy: note that each in- 
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variant polynomial F (simply called invariant) can be 
written as a sum of its homogeneous components, each 
of which is an invariant as well. One can therefore always 
find a generating set of the invariant algebra which con- 
sists of homogeneous invariants only. Furthermore, the 
set of homogeneous invariants of fixed degree is a finite- 
dimensional vector space, as one can easily verify (which 
gives the algebra of invariants the structure of a graded 
algebra). Therefore, a natural approach to our problem 
is to consider these spaces of homogeneous invariants de- 
gree by degree and to construct a basis of invariants for 
each degree. This construction will yield a generating 
(yet infinite) set of the invariant algebra. 



II. THE LOCAL CLIFFORD GROUP 



The Clifford group C\ on one 
group of unitary 2x2 matrices: 



qubit is the following 
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The order of C\ is finite and equal to 192. Up to overall 
phase factors, the Clifford group consists of all unitary 
operators which map the Pauli group to itself under con- 
jugation; here, the Pauli group Q\ (on 1 qubit) consists 
of the identity ctq and the three pauli matrices 



all having 4 possible overall phase factors equal to ±1 
or ±i. In other words, up to these overall phase factors, 
the group C\ is the normalizer of Q\ in the unitary group 
U(2). Note that these phases are not relevant in the fol- 
lowing, since we are considering the action of the Clifford 
group under conjugation as explained in the introduc- 
tion. It follows that every U £ C\ is, for our purposes, 
completely described by a permutation tt G S3, where S3 
is the symmetric group on 3 letters, and a set of three 
phases a\, 0:3 = ±1, such that 

ISo^ = OLiO^ (i = 1,2,3). 

Moreover, since <j\<ji ~ 03, one has ol\ol%olz = 1 and it is 
therefore sufficient to keep track of only two of the afs 
(say a± and a^). Another useful characterization of the 
Clifford group is obtained by considering the mapping 



a Clifford operation corresponds to an invertible linear 
transformation Q £ GL(2,¥2) (instead of a permutation 
tt G S3) and a couple of phases a\ and 0:3. It is this sec- 
ond description of Clifford operations in terms of binary 
linear transformations which is most often used in the 
literature in quantum information theory and quantum 
computing, and we will do the same. 

The local Clifford group C l n on n qubits is the n-fold 
tensor product of C\ with itself, i.e. 



C n - Ci 



) C\ (n times). 



Analogous to the case of one qubit, the group C l n can be 
most easily described by its action on the Pauli group Q n 
on n qubits, defined by 

Qn = Oi ® • • • <8> Gi (n times). 

Using the mapping the elements of Q n can be repre- 
sented as 2n-dimensional binary vectors as follows: 
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where (u,v) = (u±, . . . , u n , vt, ■ ■ ■ , v n ). As in the case 
of one single qubit, local Clifford operations map Q n to 
itself under conjugation. Therefore, n-qubit local Clifford 
operations as well can be described in terms of linear 
operations over F2. One can readily verify that, in this 
binary picture, an operator U £ C l n corresponds to an 
invertible 2n x 2n binary matrix Q of the block form 



A B 
C D 



where the n x n matrices A, B, C, D are diagonal, and a 
set of 2n phases cti = ±1, defined by 



a Qei, 



(2) 



where e; is the ith canonical basis vector in F| n , for 
every i — l,...,2n. Denoting the diagonal entries of 
A,B,C,D, respectively, by a,, 6,, c,, d,-, respectively, the 
n submatrices 



Q 



w ._ 



£ GL(2,¥ 2 ) 



correspond to the tensor factors of U. The group of all 
such Q is isomorphic to GL(2,F2) n (and S3). 
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= 011 1- 


- (1,1), 



(1) 



which establishes a homomorphism between the groups 
Qi and F|. Here, F2 is the finite field of two elements 
(0 and 1), where arithmetics are performed modulo 2. 
In this representation of Pauli matrices by pairs of bits, 



III. INVARIANT POLYNOMIALS AND 
MATRIX ALGEBRAS 

Let {pij} be a set of 2 2n variables, which are assem- 
bled in a 2™ x 2™ matrix p — (pij). Any homogeneous 
polynomial F(p) of degree r £ No can be written as a 
trace 

F(p) = Tr (A F ■ p® r ) 
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for some complex 2 nr x 2 nr matrix Ap. To see this, sim- 
ply note that the tensor product p® r contains all mono- 
mials of degree r in the entries pij. The coefficients of 
these monomials in the polynomial F are encoded in the 
entries of Ap (note, however, that the correspondence 
F <-> Ap is not one-to-one). It can easily be verified that 
F(U pU^) — F(p) for every U £ C l n if and only if there 
exists an Ap such that 



U^A F {U^ r ) 



A t 



(3) 



for every U £ C l n . Therefore, the study of invariant ho- 
mogeneous polynomials of fixed degree r is transformed 
to the study of the algebra A n , r of matrices Ap which 
satisfy ©. In this section, we will construct a linear basis 
of this algebra. This will yield a g enerating set of homo- 
geneous invariants of degree r [l9j . First we consider the 
simplest case of one single qubit, i.e. n — 1, and then we 
move to the general case of arbitrary n. 



A. One qubit 

Let r £ No and let TZ r be the averaging operator which 
maps a 2 r x 2 r matrix A to 



|Cl1 ^ 



Note that TZ r is the orthogonal projector of the space of 
2 r x 2 r matrices onto the subspace Ai, r - Therefore, a 
spanning (though in general non- minimal) set of Ai. r is 
obtained by fixing a vector space basis of 2 r x 2 r matrices 
and calculating its image under TZ r . In this context, a 
natural choice for such a basis is the set {o~r u>v \ \u,v £ F£} 
of Pauli operators on r qubits (all having an overall phase 
equal to 1). Before calculating the images TZr(&(u,v)) m 
lemma 1, we need some definitions: firstly, let the group 
GL{2,¥ 2 ) act on Ff r as follows: 

(4) 



Q £ GL(2,¥ 2 ) : {u,v) £ Ff ^ ( u ,v) £ 



where (u, v) is defined by 



Q 



for every j — l,...,r, where Uj,Vj,Uj,Vj, respectively, 
are the components of u, v, u, v. Secondly, let the binary 
vector space V r consist of all (u, v) £ F| r such that 

r 

5^(tt i ,w i ) = (0,0). 

We are now in a position to state the following lemma: 

Lemma 1 Let r £ Nq. Let (uq,vq) £ ¥ 2 r and denote by 
r the orbit of this vector under the action ^j). Then 

otherwise, 
where c is a constant. 



(«o,«o) 



Proof: Let U S C\ be an arbitrary Clifford operation. 
The action of U on the Pauli matrices is parameterized by 
coefficients aoi, oliq, olu — ±1 with aoiaioan = I and 
a linear operator Q £ GL(2, ¥2) such that U<rt at b)W = 
&abO'Q(a,b) for every (a, b) £ F^O}. Defining the integers 
n Xl n y ,n z by 

n x = \{j I (uoj,v 0j ) = (0, 1)H , 

n V = \{j I (lM)j,VQj) = (1,1)}|, 
n z = \{j I (uoj,v 0j ) = (1, 0)1| , 

the operator U® r maps o-( UOtVo s to 



T (u ,v ) 



a m a io <T (oo,«o) 



(5) 



under conjugation, where (uo,Uo) G T is the image of 
(uo,Uo) under the action Q of Q. The crucial observa- 
tion is now that the coefficient of cr( Qoi0[1 ) in (jSJ) is always 
positive (and thus equal to 1) if and only if both the 
numbers n x + n y and n z + n y are even. Note that this 
occurs if and only if n x , n y and n z are all even or all odd 
or, equivalently, if and only if (uo,^o) £ V r , as one can 
readily verify. It follows that 



(U ,V )) 



(«,u)er 



if (uq,vq) £ V r . If (uoi^o) t V r , one can easily see that 
the different terms in the sum TZr{(Xf UOjVo \) interfere such 
as to yield zero. This ends the proof. □ 
Using the result in lemma 1 , we can construct a basis of 
Ai, r - Denote by O r the set of all orbits T of the elements 
in V r (note that O r forms a partition of V r ). For every 
r £ O r , define the matrix 



At 



(6) 



By construction, the matrices Ar linearly generate the al- 
gebra Ai, r - Moreover, this set of matrices is linearly inde- 
pendent: indeed, this follows immediately from the linear 
independence of the Pauli operators cri u , v ) ■ Therefore, we 
can conclude that the Ay 's are a basis of A\, r ■ In order to 
calculate the dimension \O r \ of A\^ r , we use the Cauchy- 
Frobenius orbit-counting lemma, which states that the 
number of orbits of a finite group G acting on a set X 
is equal to the average number of fixed points, i.e., the 
number of orbits is equal to 



^Ei Fix (s)i> 

11 sec 



(7) 



where |Fix(g)| is the number of fixed points in the set 
X of the group element g. Let us therefore calculate 
the number of fixed points of an arbitrary matrix Q £ 
GL(2,¥2) acting on V r - Firstly, it is trivial that the 
identity has |V r | = 4 r_1 fixed points. Secondly, there 
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are three elements in GL(2,F2) of order two. Consider 
e.g. the matrix 



Qo = 



1 

1 



When acting on F|, this operator fixes exactly two vec- 
tors, namely (0,0) and (1,1). Therefore, when Qq acts on 
V r , the set Fix(Qo) consists of all vectors of the form 

ai(l,0,...,0; 1,0,. ..,0) + a 2 (0, 1, . . . , 0; 0,1,. ..,0) 

+ ... + a r (0,0,...l; 0,0,..., 1), (8) 

where on £ {0, 1} for every i = 1, . . . , r and where ex- 
actly an even number of a^s are nonzero. Therefore, 
the cardinality of Fix(Q ) is equal to the number of even 
subsets of {1, ... , r}, i.e. |Fix(Q )| = 2 r ~ 1 . Note that an 
analogous argument holds for the other two matrices of 
order two. Finally, there are two elements in GL(2,¥2) 
of order three, which fix only the zero vector. Gathering 
these results in the formula Q, we find that the number 
\O r \ of orbits is equal to 



1 



:(4 r 



3-2 r 



2). 



We have proven: 



Theorem 1 Let r £ Nq. The set {Ar}reo r is a vector 
space basis of the algebra Ai >r . The dimension \O r \ of 
A\ t r is equal to 



-(2 



2r ~ 3 + 3 • 2 r ~ 2 



1). 



(9) 



Thus, we have obtained the desired result of construct- 
ing a basis of matrices of the algebra Ai <r . It will be 
useful to have an explicit parameterization of the orbits 
r £ O r . Such a parameterization could e.g. be used to 
enumerate all the matrices A? for a given degree. Also 
when we will move from the matrix algebra Ai. r to the 
polynomials Tr (Ar-p® r ) in section 4, a more operational 
description of the A-p's will turn out to be very useful. 
To this end, for each (u, v) £ F 2 , 7 ", define the sets 

r) (u,v) = {j\ (uj,Vj) = (0,0)}, 
Vx(u,v) = {j\ (Uj,Vj) = (0,1)}, 
T)y(u,v) = {j\ (Uj,Vj) = (1,1)}, 

Vz(u,v) = {j\ (u^Vj) = (1,0)}. 

Then the following characterization is easily verified: two 
vectors (u,v), (u',v') £ W 2 , 7 " belong to the same orbit if 
and only if 

(a) rjo(u,v) = r) (u',v') and 

(b) there exists a permutation n of {x,y, z} such that 
n x (u',v') = )7 ff ( x )(u,u), r) y (u',v') = r) v ( y )(u,v), and 
r) z (u',v') = r)^g)(u,v). 

This implies that any orbit T of the action (JIJ can com- 
pletely be described by 



(a') a set rjo(T) C {1, . . . , r} and 

(b') a partition V(T) = V2 , m } of {1, . . . , r} \ ry (r) 
into three (possibly empty) subsets, 

such that (it, v) £ T if and only if r)a(u,v) = r/o(T) and 
{r) x (u, v), r] y (u, v),rj z (u, v)} = V(T). Moreover, T £ O r if 
and only if the numbers |ryi|, 1 772 1 , 1 773 1 are either all even 
or all odd (cfr. proof of lemma 1). Let us illustrate this 
characterization with two simple examples: 

• r = 1: there is one orbit in 0\, namely Tq = 
{(0,0)} £ 0\. This orbit is characterized by 
ry (r ) = {1} and V(T ) = {0,0,0}. 

• r = 2: there are two orbits in O2, namely T — 
{(0,0; 0,0)} and 

r' = {(0,0; 1,1), (1,1; 0,0), (1,1; 1,1)} 

= {(«,«) € Ft|(«i,«i) = (u 2 ,v 2 ) + (0,0)}. 

The orbits T and T' are described by 

?/0 (r) = {l,2}, V(T) = {0,0,0} 

and 

77o(r') = 0, P(r') = {{i,2},0,0}. 

B. Multiple qubits 

For arbitrary n, the result in theorem 1 can immedi- 
ately be used to construct a basis of A n , r - To see this, 
let us first consider the algebra of 2 nr x 2 nr matrices A 
which satisfy 



Uf 



U® r A{Uf 



u 



5r\f _ 



.4. 



for every U±, ...,[/„ £ C\. It is straightforward to show 
that this algebra is the n-fold tensor product of A\, T with 
itself. Therefore, a basis of this algebra is given by the 
matrices Ar 1 (8 • • • <8> Ar n , where r, ranges over all orbits 
in O r , for every i = 1, . . . , n. In order to obtain a basis 
of A n ,r, one simply has to conjugate this basis with the 
permutation matrix P, defined by 

P \i\l ■ ■ ■ Hr] «21 ■ ■ • «2r; • • • 5 «nl • • ■ w) 

= |«11 • • • inV, il2 ■ ■ ■ in2\ ■ ■ ■ \hr ■ ■ ■ w), (10) 



where i a b £ {0,1} and |*iix • - * ) are t ne standard ba- 
Indeed, the matrix P performs 



sis vectors in 



12"*" 



the appropriate permutation of tensor factors, mapping 
Uf r ® ■ ■ ■ <g> U® r to{U x ®---® U n )® r under conjugation. 
This leads to the following result: 

Theorem 2 Let r £ N. For every n-tuple 7 = 
(ri, . . . , r„) of orbits Ti £ O r , define the matrix 



-4, 



PA 



A v P T . 



ill) 



Then the set {Ay} 7 forms a vector space basis of A n . r 
The dimension of A n ,r is equal to \O r \ n . 
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Following the discussion at the end of section 3.1., the 
matrices Ay can be described in an alternative way than 
lfTT)l . using the description of orbits T £ O r by couples 
(770 (r), V(T)). Defining the support of a vector w £ F|™ 
to be the set 

supp(w) ={ie{l,...,n} I (Wi,w n+i ) ^ (0, 0)}, (12) 
one obtains: 

Theorem 3 Let 7 = (Ti, . . . , r„) be an n-tuple of orbits 
Tt £ O r - For every j,k £ {1, . . . , r}, j < fe, define the 
sets u/(fi and 67/ 



/J) 



{ie{l,...,n} I j £ r) (Ti)} 
w (i fe ) = {i e {1, . . . , n} I j, fc € f?o(rj) or j and fc 

belong to the same subset ofV(Ti)}. (13) 

Then A~ i — ^2 o~ w (i) (8- • •<8o~. u ,(r) , where the sum runs over 
all ordered r-tuples (w^\ . . . ,w^) £ (F 2 , 11 )* 1, satisfying 



supplw^) = u;V' 

SUpp(lU^ + W^) = LO^ k \ 



(14) 
(15) 



for every j,k £ {1, . . . ,r}, j < k, where uj^\ ui^ k ^ de- 
note the complements of the sets 

U U) > U W m{l,...,n}. 
Proof: By definition, A 7 is equal to 

^2 a wW ® • • • ® cr ro (r) , 

where the sum runs over all ordered r-tuples 
. .,ioW) € (F^ l ) xr such that 



(w. 



(i) 



W (i) 



n+i/ c 



(16) 



for every i = 1, . . . , n. The proof of the theorem then fol- 
lows immediately from the characterization of the orbits 
Ti by the couples (7/0(1^), V{Ti)), for every i = 1, ... , n. 
□ 

Example 1. Let us consider this result for the case 
of smallest nontrivial degree, i.e. r = 2. Let 7^) = 
(ri, . . . , F„) be an n-tuple of orbits L G C2- Recall that 
02 contains exactly two orbits Y and T', as defined in 
the last paragraph of section 3.1. Let 10 be the subset 
of {1, . . . ,n} which consists of all i such that F^ = F. 
Following the definitions stated in theorem 3, we have 



UJ 



i) 



: UJ 



/ 2 ) and a/ 12 ' = {1, . . . ,n}. Consequently 



Ay(2) 



uj£F| n , supp(u;)— a) 



This shows that the matrices A ( 2 ) are parameterized by 
the subsets w of {l,...,n} in a one-to-one correspon- 
dence. 

While the result in theorem 3 is in fact no more than 
a reformulation of (|11|) . it is interesting in that it relates 
the matrices Ay (and thus the corresponding invariant 
polynomials Tr (Ay • p® r ) as well) to the notion of the 
support of a binary vector, which is of central importance 
in quantum coding theory. Note that the definition l|12|) 
of support is indeed the same as is used in the theory of 
quantum codes. 



IV. BASES OF INVARIANTS 



It follows from theorem 2 that the polynomials 



plM :=Tr(Ay^), 



(17) 



in the variables pij (i,j = 0, . . . , 2™ — 1) linearly gener- 
ate the space of homogeneous invariants of C l n of degree 
r. However, different Ay's may correspond to the same 
polynomial and therefore linear dependencies within the 
set of the polynomials H17|) can exist in general. We now 
set out to pinpoint a basis of polynomials for each degree 
r. As in the preceding section, we start by considering 
the simplest case of one qubit and then move to the gen- 
eral case. 



A. One qubit 

Let p = (pij), where i,j — 0,1, be a matrix of vari- 
ables. Fix an orbit T £ O r with r/o(r) = 770 and 
V(T) = {771,772,773}. It will be convenient to introduce 
the linear forms Xij(p) := Tr(/5<xy), where i.j = 0, 1, or 
more explicitly: 



xoo(p) = Poo + Pn 

xoi(p) = P01 + Pw 

%w(p) = Poo - pn 

xn(p) = Kpoi ~ Pw)- 



(18) 



Conversely, the pij 's can be written as linear forms in the 
variables x — (xij) as follows: 



I 1 



»,i=o 

We will consider Tr (A^ • p(x)® r ) to be a polynomial in 
the variables x. This yields 



Tr (Ar • p(x)® r ) = 1 £ x UlVl 

= L E 

or 



(u,»)er 



2r 00 01 

(u,v)ET 



ns(u,v) n^(u,v) n z (u,v) n y (u,v) 



Xr\n Xn-i Xm X 



'11 1 



(19) 



where we have used the definitions no(u,v) = \f]o(u, v)\ 
etc.. Note that 



and 



no(u,v) = \rjo\ 



{n x (u,v),n y (u,v),n z (u,v)} = {|ni|, |77 2 |, |7? 3 |} 



for every (u, v) £ T. It readily follows that l|19l) is equal 
to 



J>7ol V 



Jl«(l)U?>WlJl.(3)l 

x Q1 x 1Q x 11 



(20) 
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up to a normalization factor. Expression 1|20|) shows that 
the polynomial Tr (A-p ■ p® r ) only depends on the number 
| »7o | an d the set {|ni|, 1 772 1 , 1 773 1 } . In other words, if T and 
r' are two orbits such that 

|7»(r)| = Mr')| 

and 

{Mr)|, \ m (r)\, \ m (T)\} = {Mr')|, |772(r')|, \m(r')\}, 

then (and only then) the polynomials Tr (Ar ■ p® r ) and 
Tr (Ar> ■ p® r ) coincide. This equivalence relation on 
O r leads to the following definition: for each 4-tuple 
A = (Ao, Ax, A2, A3) of non-negative integers A^ such that 
Ai, A2, A3 are either all even or all odd, A0 + A1 + A2 + A3 = 
r and Ai > A2 > A3, we define an invariant p^ of C\ of 
degree r as follows: 

Pr = ^00 x oi ( ' x io ( ' x n ( ' ■ (21) 

Recall that p* is to be regarded as a polynomial in the 
variables p via H18|) . By construction, the set of all these 
polynomials generates the space of invariants of degree 
r. What is more, the p^ 's are linearly independent. This 
immediately follows from the fact that each monomial 
in the variables Xij occurs in exactly one polynomial p* 
and that the polynomials Xij(p) are algebraically inde- 
pendent. We have therefore proven: 

Theorem 4 The polynomials p^: form a basis of the vec- 
tor space of homogeneous invariants of C\ of degree r. 

B. Multiple qubits 

The construction of bases of invariants for arbitrary n 
will be a generalization of the one qubit case. Starting 
from a 2™ x 2™ matrix p of variables, we again perform 
a change of variables, defining x w = x w (p) — Tr(p • a w ), 
for every w G F| n . Analogous to the one qubit case, the 
converse relation reads p(x) = i Ylw x w°~w Note that 
the polynomials {x w (p)} are algebraically independent; 
this follows from the fact that the variables x and the 
variables p are related by an invertible linear transfor- 
mation. Now, letting 7 = (T\, . . . , T n ) be an n-tuple of 
orbits G O r , the invariant p^ r , regarded as a polyno- 
mial in the variables x, is equal to 

x wW ...x w( r) (22) 

(wC-t ,...,!!;(»•)) £7 

up to a normalization. Here, (u/ 1 ), . . . , w^) G 7 is a 
shorthand notation to express that (ai' 1 ',...,/)) is an 
r-tuple of vectors w^> G F2™ satisfying 

(wP,...,w < f\w^_ il ... 1 w^eT l , (23) 

for every i = 1, . . . , n. As in the case of one single qubit, 
the correspondence between the polynomial p^ n r and the 



matrix A y is non-unique. Indeed, suppose that p G S r 
is an arbitrary permutation and define the n-tuple 7^ = 
(It,...,r£) such that 

jevaW) iff /^COe^r*) (24) 

for every j G {l,...,r} and a G {0,1,2,3}. Equiv- 
alently, one has (w^, . . . , w^) G 7 M if and only if 
(w < -^\...,w { ^>) G7. Then 

Pn,r = Pn,ri (25) 

which immediately follows from (|22|l . Conversely, if 7 
and 7' are two n-tuples of orbits such that p 1 nr = p^ r , 
then there exists a permutation p G S r such that 7' = 
7 M , as one can easily verify. We now claim that a basis 
{pZ] r :Pn 2 ,r: ■ ■ ■ } of the space of invariants of is obtained 
by fixing a set {71 , 72 , . . . } of n-tuples of orbits such that 

(i) The polynomials p^ r are pairwise different 

(ii) For every n-tuple 7 of orbits, p'^ r = p"J* T for some 
i = l,2,.... 

The claim is proven as follows: firstly, it follows from 
the construction of the invariants p% r and item (ii) that 
the polynomials p^ r generate the space of homogeneous 
invariants of degree r. Secondly, the linear independence 
of the p2v' s follows from (i). For, suppose there exist 
complex coefficients a,, not all equal to zero, such that 

= 0. (26) 

i 

As each monomial n^=i x wtii 1 where 10 W G F2' 1 , occurs 
in exactly one invariant p^ r , this yields a nontrivial linear 
combination of these monomials adding up to zero, which 
is a contradiction; indeed, the monomials Oj=i x w<J) are 
linearly independent, as the polynomials {x w (p)} are al- 
gebraically independent. 

We now set out to construct a set of invariants which 
satisfies (i)-(ii). According to the discussion above, there 
is an equivalence relation ~ on the set 0™ of n-tuples 
of orbits, such that 7 ~ 7' if and only if there exists a 
permutation p, € S r such that 7' = 7^. A set of invari- 
ants which satisfies the desired conditions is obtained by 
choosing any set {71 , 72 , . . . } of orbits such that every 
equivalence class is represented by exactly one n-tuple 

Recall that an n-tuple 7 = (Ti, . . . ,T n ) G 0™ is de- 
scribed by n couples (ryo(r,), V(Ti)), where T]o(Ti) C 
{1, . . . ,r} and V{Ti) is a partition of {1, . . . ,r} \ f]o(^i) 
into three subsets. While such a system of n couples 
compactly describes 7, it will be useful to represent 7 in 
a different way, which contains some redundant informa- 
tion but has the advantage of being more transparent: 
we describe 7 by an n x r matrix M with entries in the 
set {0, 1, 2, 3}, satisfying 

Afy=0 iff j G? ? o(lU 

7^ Mij = Mik iff j and k belong to the same subset 
in the partition V(Ti), (27) 
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for every i = 1, . . . , n and j, k = 1, . . . , r. It is clear that 
this description exhibits some degeneracy, as any permu- 
tation of {1,2,3} in any row of M yields a (generally) 
different matrix which also satisfies l|27|). However, the 
equivalence relation ~ is translated into a simple kind 
of equivalence transformation of matrices. Indeed, two 
n-tuples 7,7' € 0™, described by n x r matrices M and 
M', respectively, belong to the same equivalence class 
of the relation ~ if and only if M' is equal to M mod- 
ulo a permutation [i S S r of its columns and n row-wise 
permutations iii of {1, 2, 3}, and we write M ~ M '. 

Seeing that we are looking for suitable representatives 
of each equivalence class, it is appropriate to look for 
normal forms of the matrices M under the above action 
of the permutations (i and 7Tj. There is in fact a lot of 
freedom to define sensible normal forms. One possible 
definition is stated below in definition 4. First we need 
some preliminary definitions: 

Definition 1 Let d G Nq. Let u — (ui, 1*2, • • • , Ud), 
v = (t>i, x>2, ■ ■ ■ , Vd) be two d- dimensional vectors with 
nonnegative integer components. A lexicographical order- 
ing relation <i ex is defined as follows: u <i ex v if it = v 
or if there exists j (1 < j < d) such that it, = Mj if i < j 
and Uj < Vj . 

Definition 2 Let u be a d- dimensional vector with en- 
tries in {0,1,2,3}. For every a G {0,1,2,3}, define 
Va(u) = {j G {1, . . . , d} I Uj = a}. 

Definition 3 Let M be an n x r matrix with entries 
in the set {0,1,2,3}. Let Mf denote the ith row of 
AI . Let m = (mi, . . . , mj D ) be an iq- dimensional vec- 
tor with entries in {0,1,2,3}, where io < n. Then the 
set r) m (M) C {1, . . . r} is defined as follows: 



r} m (M) = C\i< io ri mi { M i)- 



(28) 



For every a G {1,2,3}, the vector uf\.\{M) with 

components u^ +1 (M) m , where m ranges over all ig- 
dimensional vectors with components in {0, 1, 2, 3}, is de- 
fined by 

«tli(% = \{j G Vm(M)\M i0+hj = a}\ (29) 

(the indices m of the components of u(^ +1 {M) are ordered 
according to the lexicographical ordering relation.) 

Definition 4 Let M be an n x r matrix with entries in 
the set {0, 1, 2, 3}. Then M is in normal form if it satis- 
fies the following conditions: 

(i) The columns Kj of M are ordered non- deer easingly, 



i.e. Ki <i 



^ l.P.T. K r 



(ii) \m{Ml)\ < \V2{M?)\ < \m{M?)\ and for every 
i = 2,...,n, 



(iii) For every i = 1, ...,n the three numbers 
\r)i(MT), \i] 2 (Mf )\, \tj 3 (MT)\ are either all even 
or all odd. 

Example 2. The following 3 x 11 array is in normal 
form: 



111 1 22 33 

12 111 2 33 22 

1 2 3 012 3 03 12 



(31) 



Indeed, conditions (i) and and (iii) are easily checked, as 
well as the first part of condition (ii) . As for the second 
part of (ii), let us calculate the vectors 

ui a) (M) = ((4 a) ) , (ui%, (u { 2 %, (ui%) (32) 
and (M) is equal to 

((4 a) ) 00 , (4 0) )oi> (4 a) ) 02 , (4 a) )o3, (4 a) )io, (4 Q) )n, • ■ ■) 

Using definition (|29|l . we find 



,(!) 



and 



(1,3,*,*), 4 = (1,1,*,*), u 



(3) 



(0,0,*,*) 



4 X) = 
4 2) = 

,(3) _ 



(1,0,0,*,...) 
(0,1,0,*,...) 
(0,0,1,*,...), 



(33) 



uf\M) < lex 4 2) (M) <ie X uf\M). 



(30) 



where the entries denoted with * are (in this example) 
irrelevant to order the vectors lexicographically, and con- 
dition (ii) follows. o 
One can easily verify that each equivalence class con- 
tains exactly one normal form. Note that, given an n x r 
normal form M, one recovers the corresponding tuple 
7m = (Ti, . . . , T„) G 0" as follows: 

ViTi) = { m (M^), m (M[), m (Ml)} . (34) 

For instance, the tuple 7 corresponding to the normal 
form in example 2 is defined by: 

77o(ri) - {1,2,3}, P(r 1 )-{{4,5,6,7},{8,9},{10,11}}, 
J7o(r 2 ) - {1}, V(T 2 ) ={{2, 4, 5, 6}, {3, 7, 10, 11}, {8, 9}}, 
77o(r 3 ) - {4,8}, P(r 3 )-{{1,5,10},{2,6,11},{3,7,9}}. 

We have proven our main result: 

Theorem 5 For every nxr normal form M , denote the 
corresponding n-tuple of orbits by 7m . Then the set of all 
invariants p^ M r forms a basis of the space of homogeneous 
invariants of C l n of degree r. 
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Thus, we have obtained our initial objective of construct- 
ing for every n and for every r a basis of the space of 
invariants of C l n of degree r. Note that for the case n = 1 
we indeed recover the result obtained in the previous sec- 
tion. 

It is interesting to investigate the behavior of the di- 
mensions d n ^ r of these spaces for large n and r. Lower 
and upper bounds for d n ^ r are the following: 

Lemma 2 Let n, r e No- Then 



we find that 



6™r! 



(4 



r-1 



3 • 2 r_1 + 2) n < d n ,r < 



4"_ 1 
r 



Proof: Let M. n y.r denote the set of all n x r matrices 
M with entries in the set {0, 1, 2, 3}, such that for every 
i — 1, . . . , n the three numbers 



MM?), \m(M?)\, MM?)] 



(35) 



are either all even or all odd. Recall that d„ lT - is equal 
to the number of orbits of the group S r X S3 1 acting on 
this set as defined above. Using the Cauchy-Frobenius 
lemma, the number of orbits is equal to 



— Fix(/Z,7R), 



(36) 



where Fix(/x,7Ti) denotes the number of fixed points in 
■M nX r °f the element (//,7Tj) = (fx, 7Ti, . . . , 7T n ), where 
/i G <SV and 7Tj € S3. Firstly, note that restricting the 
sum to all group elements were /i is equal to the identity 
yields the desired lower bound, using a highly similar 
argument to the one used to calculate \O r \ n above. In 
order to obtain the upper bound, we will calculate the 
number N n>r of orbits of the group S r acting on the set 
of all n x r matrices with entries in the set {0,1,2,3} 
by permuting columns. Note that this number is indeed 
an upper bound for d n ^ r . The Cauchy-Frobenius lemma 
yields 



h E (4": 



c(/i) 



(37) 



where c(/i) denotes the number of cycles in the permuta- 
tion jj,. Consequently 



(38) 



k=Q 



where t{r, k) is defined as the number of permutations in 
S r which have exactly k cycles. Note that this number 
is related to the Stirling number s(r, k) of the first kind 
by the relation t(r,k) = (—l) r+k s(r,k) [16j. Using the 
identity [HI 



fc=0 



s{r,k)x k = (-l) r r! 



x 
r 



(39) 



r + 4" - 1 
r 



(40) 



which completes the proof. □ 
While these bounds are in fact quite rough, they are 
sufficient to gain qualitative insight into the limit behav- 
ior of the dimensions d n . r when n or r are large. Let us 
first examine lim r ^oo d n ,r for fixed n. Denote A = 4" — 1. 
Then, using the Stirling approximation ln(a!) ~ a In a— a, 
the upper bound reads 

ln^ r + A ^ = ln(r + A) ! - In r! — In A! 

w (r + A) ln(r + A) - r In r - In A! - A 

= ln(l + -) r + Aln(r + A) - In A! - A 
r 

w Aln(r + A) - In A!, (41) 

where in the last line we have used (1 + ^) r w exp(A) 
when r is large. Finally, we obtain 



d n ,r < ^(r + A) ; 



(42) 



We have proven: 

Theorem 6 For every fixed n £ No, the dimension d n ^ r 
tends polynomially in r to infinity . In other words, for 
every n there exists a polynomial p n (r) in r such that 
d n ,r = O (jp n (r)). 

Note that a similar result does not hold for limn—xx, d n , r 
for fixed r. Indeed, the lower bound in lemma 2 shows 
that 



which is nonpolynomial in n if r > 2. 



(43) 



V. INVARIANTS OF DEGREES 1, 2 AND 3 

In this section we investigate the invariants of C l n of 
low degrees in more detail. In particular, we will show 
the following result: 

Theorem 7 Every invariant of C l n of degree 1, 2 or 3 is 
an invariant ofU(2)® n (which also acts by conjugation) 
and vice versa. 

One of the implications in the theorem is trivial. Indeed, 
every invariant of U(2)® n is an invariant of C n , as the 
latter is a subgroup of the former. Let us now prove the 
reverse implication. 

Let p be a 2™ x 2™ matrix of variables. Firstly, it follows 
from theorems 1 and 2 that C l n has only one invariant of 
degree 1, namely Tr(p), which is trivially an invariant of 
U{2)® n . 

In order to examine the invariants of degrees 2 and 3, 
it will be convenient to introduce the following functions: 
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Definition 5 Let lj C {1 , . . . , n}. Define the functions 
:F|"^C by ' 

Sw(w) — 1 if supp(w) = lo and 5 u) {w) — otherwise 
£uj(w) = 1 if supp(w) C lo and e UJ (w) = otherwise. 
It is straightforward to show the following relations 

e u = E S u i 

Lj' C.LJ 

S u = (-I)'"' E (-l) 1 "' 1 e<y, (44) 

the first of which is trivial and the second of which can 
easily be verified by substitution in the first one. 

Now, regarding r = 2, using example 1 we find that 
the polynomials 

PUP) = E Tr (a w <S> a w ■ p® 2 ) 

E Tr{(a w -pf} (45) 

supp(u;)— u) 

where lo ranges over all 2™ subsets of {1, . . . ,n}, form 
a generating set of the space of invariants of degree 2. 
Moreover, using the techniques of the previous section, 
one can easily show that the p u 's are linearly independent 
and therefore the dimension of this space is 2 n . Interest- 
ing variants of 1 (45( 1 are the polynomials 

qu>(p) = E Tr {(cr^ • p) 2 } 

tu£F| n , supp(iu)Ca; 

= Tr {(Tr* p) 2 }, (46) 

where the operation Tr^ denotes the partial trace over all 
qubits outside the set lo. The polynomials q u are mani- 
festly invariant under the entire local unitary group. In 
fact, it is well known that these polynomials are gener- 
ators of the space of invariants of U(2)® n of degree two 
[l7j . Moreover, one has the relations 

1" = E Plo ' 

Lo'C-LO 

p» = (-i) H E ( 4? ) 

which follow immediately from (|44|l . In particular, the 
second expression in 1)47(1 shows that every polynomial 
Puj is an invariant of [/(2)® n , implying that the sets {p^} 
and {qui} span the same space, which yields the desired 
result for theorem 6 for r = 2. Furthermore, it follows 
from (|47|l that polynomials q w are a basis as well, being 
a generating set of cardinality 2" in a 2"-dimensional 
space. 

A similar result can be proven for the invariants of 
degree 3. Theorem 2 shows that the space of invariants 
of C l n of degree 3 is spanned by all polynomials 

Pn,3 = E Tl ' ^ w{1) ® a w^) ®°~w(.3) P m ), 

(w( 1 ),to< 2 ),i«( 3 ))e7 



where 7 ranges over all elements in O3 . Note that, for 
every 7 £ O3 , one has it/ 1 ' + 1// 2 ' + = whenever 
(it/ 1 ',-!// 2 ',-!// 3 ') S 7, by definition of O3. Using the de- 
scription of 7 by sets o;W and w^'' introduced in theorem 
3, it follows that 

<3 = E Tr ® CT ™ <2> ® ff wM+v,W P® 3 ), (48) 

where the sum runs over all couples (w^\w^) G 
(F 2 .")* 2 such that 

SUpp(w^ 1 ^) = LUl, SUpp(u/ 2 ') = Ll>2 

supp(w (1) + w (2) ) = W12 , (49) 

for some uj\, 002, 0012 Q {1, . . . ,n}. Using 1 (11)1 . a straight- 
forward calculation shows that p^ 3 is, up to an overall 
sign, equal to 

E (-I)^l + '^l+K 2 l Tr {(t M/9) (Tr _, p) (Tr^p)}, 

(50) 

where the sum runs over all C wj, U2 C W2 and 
w 12 C W12. As the summands in 1)50(1 are manifestly in- 
variant under the action of U(2)® n , the polynomial p„ 3 
is an invariant of the local unitary group and the proof 
of theorem 7 is completed. 

VI. CONCLUSION 

We have performed a systematic study of the invariant 
algebra of the local Clifford group C l n , using the descrip- 
tion of this group in terms of binary arithmetic. Our 
approach was to consider the spaces of homogeneous in- 
variants degree per degree and to construct bases of these 
spaces for each degree r. In order to study these spaces of 
homogeneous invariants, we transformed the problem to 
the study of certain algebras A n . r of matrices, such that 
every matrix in an algebra A n . r corresponds to an invari- 
ant polynomial of degree r. We then constructed bases 
{Ay} 7(E 0;.* of these algebras, which yielded generating, 
though linearly dependent, sets {p^ r } 7 of homogeneous 
invariants. We subsequently showed how a basis of in- 
variants could be pinpointed amongst these polynomials 
for each degree r, which was the main result of this paper. 

As stated in the introduction, we believe that these 
results are relevant in a number of fields in quantum in- 
formation theory, with in particular, the classification of 
binary quantum codes. In forthcoming work we will ap- 
ply the present results to this problem. 
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